First order phase transition in the quark matter 
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We investigate chiral phase transition of the first order in the quark matter. Using the Nambu- 
Jona-Lasinio model, an equation of state of the quark matter which is similar to the van der Waals' 
one is obtained. Moreover the specific heat and the compressibility are calculated. It is shown that 
they are enhanced in the symmetry broken phase, in particular diverge near the tricritical point. 
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I. INTRODUCTION 

Thermodynamics of the quark matter have been dis- 
cussed recently by many authors One of the rea- 
sons comes from the possibility of the quark-gluon plasma 
(QGP), which lies in the chirally symmetric phase. It 
is expected that the QGP may be realized in high en- 
ergy heavy-ion collisions at CERN Large Hadron Collider 
(LHC) and at the BNL Relativistic Heavy Ion Collider 
(RHIC). On the other hand the usual hadronic matter 
is in the chiral symmetry broken phase. Therefore it is 
of fundamental interest in the study of the chiral phase 
transition in the quark matter. 

The chiral dynamics of the quark matter has been stud- 
ied mainly in the framework of the Nambu-Jona-Lasinio 
(NJL) model or QCD-like theory According to 

them, the phase transition is second order in the three- 
flavor quark matter. On the other hand, it is second 
order at low density and first order at high density in 
the two-flavor case. The critical point between them is 
called a tricritical point 0]-^5|. This behavior of the 
phase diagram was first shown by the NJL model |j| and 
then confirmed by the QCD-like model ^ is very 

interesting for us to have the first order transition. The 
present author has calculated the latent heat appeared 
in the transition and pointed out that it may be used as 
a signal of the chiral transition |16| . 

The classical liquid-gas phase transition is the typical 
first order one. There exists large fluctuation near the 
tricritical point. This phenomenon is known as "opales- 
cence" first observed by T.Andrews in the nineteenth cen- 
tury 01 • ^ is generally expected that there is also large 
fluctuation in the quark matter and some physical quan- 
tities are enhanced near the point. These may be other 
signals of QGP. It is the purpose of this paper to calculate 
such quantities, the specific heat and the compressibility 
and investigate their behaviors near the critical temper- 
ature, in particular the tricritical point. In the course of 
this work we will obtain an equation of state of the quark 
matter, which is similar to the van der Waals' one. 

The outline of the present paper is as follows. In the 
next section, we introduce the NJL model and calculate 
thermodynamic potential with the use of the mean field 



approximation. In Sec. Ill, the specific heat and the com- 
pressibility are derived from the potential. The numer- 
ical calculation of the thermodynamic quantities is car- 
ried out in Sec. IV. Section V is devoted to discussions 
and summary with the help of the Landau theory of the 
phase transition. 



II. THERMODYNAMIC POTENTIAL 

In this section we will derive the thermodynamic po- 
tential of the quark matter in the frame of the NJL model 
|18|| . Our quark matter is supposed that the number of 



flavors is N 



f 



2 and that of colors is N r 



3. It is 



known that there is the chiral transition of the first or- 
der and the tricritical point exists in this case The 
Lagrangian density of the NJL model is defined by 



C = ip(id — m)ip + g{(tpip) 2 + (ipi^Tijj) 



(1) 



where m is a current u- and d-quark mass and neglected 
in this paper for the convenience (chiral limit). The g in 
this equation denotes a coupling constant of the quark- 
quark interaction and r is the Pauli matrix for flavors. 
This Lagrangian leads to the Hamiltonian density, 



Tt = ■ V — m)ip — g{(ipijj) 2 + (^i^Tip) 2 



(2) 



Here we introduce two mean fields: a± = (ipip) and 
(72 = (tpj ip). The expectation values represent thermo- 
dynamic averages and determined from the variational 
principle later. Then we obtain the linearized Hamilto- 
nian density, 



H mf = V>(-t7 • V + M + -(U 27 ) V + Ga\ - 



2N n 



(3) 



where G = (4:N C + l)g/AN c is the renormalizcd coupling 
constant and M is the effective quark mass defined by 



M = m-2G(^). 



(4) 



Now let us introduce the thermodynamic potential 
fi(V, T, n) of our system by {0 = 1/kT) 



n 



-fcTlogTrexp [-0(H mf - fiN)] = Vco(T, /i), (5) 
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where V and T represent the volume and the temperature 
of our system respectively and y, is the chemical poten- 
tial. The mean field Hamiltonian describes an assembly 



2 



of free quasi-particles so that we can easily calculate the 
thermodynamic potential. Substituting Eq.(3) into this 
equation, we get 

cj = E a - -^T \ p 2 [log(l + exp -0(E - p r )) 
77 Jo 

+ log(l + exp —0(E + fi r ))] dp, (6) 

where a cut-off constant A is introduced as usual and the 
first term En is given by 

m a * ~ — L p Edp - (7) 



En — Go - ? — — u 2 



The quasi-particle energy and the renormalized chemical 
potential are defined by E = \/p 2 + M 2 and p r = M — 
(g/N c )<T2 respectively. This potential u plays a central 
role in the following discussion. 

Lastly we must determine the mean fields <j\ and oi. 
We take the variational principle of the thermodynamic 
potential: du/dai — du/da2 — 0. These equations lead 
to selfconsistency ones, 



_ MN c N f rA 



Jo t 1 ~ n (P>Mr) ~ m(p,p r )]dp, 



(72 = ^r 1 J A [n(p, Mr) - m(p, /v)] dp. 



(8) 



The first equation is nothing but the gap equation in the 
BCS theory and the second determines the renormalized 
chemical potential. 



III. SPECIFIC HEAT AND COMPRESSIBILITY 

Now let us calculate thermodynamic quantities of the 
quark matter in order to study the critical behavior near 
the chiral phase transition. The specific heat and the 
compressibility reflect the large fluctuation, because they 
are proportional to the fluctuations of the entropy and 
the density respectively. The thermodynamic quantities 
are classified into two groups: extensive state quantities 
and intensive ones. As for the extensive state quantities, 
we use those per unit volume in this paper. For example, 
we consider rather the number density than the volume 
itself and discuss the thermodynamic potential per unit 
volume and so on. 

The differential of the thermodynamic potential reads 

to 



duj = —SdT — pdp, 



(9) 



where S and p denote the entropy density and the num- 
ber density respectively. These quantities have disconti- 
nuities at the critical temperature of the transition |16| . 
From Eq.(9), they are represented by 



S 



N c N f 



p 2 [n log(n) + (1 — n) log(l — n) 



+m log(m) + (1 — m) log(l — to)] dp, 



N c N f 



p 2 [n — m] dp. 



(10) 
(11) 



This expression for the density is nothing but the mean 
field (72, which is a natural result from the definition. At 
this stage it should be noticed that we can obtain an 
equation of state for the quark matter. The pressure of 
the quark matter is given by p = —lu(T, p). On the other 
hand Eq.(ll) shows that p is a function of T and p. From 
these equations we get the equation of state: p = F(T, p). 

Next let us consider the second derivatives of the ther- 
modynamic potential, the specific heat and the compress- 
ibility. These quantities reflecting the fluctuation are 
supposed to be sensitive to the phase transition. They 
are defined by 



C 



T 



dS_ 
df 



lfd£ 

P \8py T 



(12) 
(13) 



However this expression for k is not useful for our cal- 
culation. It should be transformed into more tractable 
form. Using the formula of the Jacobian, it is rewritten 
to 



dp 
dp 



d(p,T) _ 3(£'t) _ Gm) t _l( Up 
d{p,T) 9 (p-T) 



9(M,T) 



p \dp 



with the use of p = —dui/dp = dp/dp. From this trans- 
formation, the compressibility turns out 



dp- 1 
dp 



(14) 



We will calculate these equations numerically and discuss 
the chiral phase transition in the next section. 



IV. NUMERICAL RESULTS 

In order to calculate the thermodynamic quantities 
considered in the previous section, the first task to do 
is to solve the self-consistency equations (8) with given 
values p and T. The parameters of the NJL model 
are the same as those used in Ref.Q: g — 5.074 x 
10~ 6 MeV -2 , A = 631MeV,TO = 0. As shown in the pa- 
per, the qualitative features of the phase transition do not 
depend on them. Then we have two solutions: M = 
and M = Mi > 0. The former corresponds to the chiral 
symmetric phase and the latter to the broken phase. An 
example with p — 0.31GeV (first order transition) is il- 
lustrated as a function of T in Fig.l. The third solution 
appears in the range Ti < T < T%. It is, however, always 
unstable so that we do not consider it. The stable state 
is realized by the solution with lower thermodynamic po- 
tential: The other solution corresponds to a metastable 
state. When the two solutions give the same value for the 
thermodynamic potential (degenerate) , the phase transi- 
tion occurs at the critical temperature T — T c . Thus we 
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FIG. 1: The effective mass M as a function of the temperature 
at fi = 0.31GeV. The critical temperatures T c , T\ and T2 are 
defined in the text. 



FIG. 3: The equation of state of the quark matter at T 
30MeV. 
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FIG. 2: The critical temperature T c as a function of the chem- 
ical potential fj,. The tricritical temperature is denoted by 
Ttcp • 



have a phase diagram shown in Fig. 2 where the critical 
temperatu re T f is drawn as a function of the chemical po- 
tential fli^-plf . The broken phase (hadron) is realized in 
the lower temperature (T < T c ) and the symmetric phase 
(QGP) in the higher one (T > T c ). The tricritical point 
is denoted by T tcp . It should be noted that the solution 
Mi with the symmetry breaking is metastable even in 
the symmetric phase above T c . This situation is known 
as supercooling. This metastable state exists in the range 
T c < T < T\ as shown by the broken line in Fig.l. On the 
other hand, the solution M — is also metastable in the 
broken phase (T2 < T <T C ). Hence the domain enclosed 
by the two dotted lines (T2 < T < T\) has metastable 
states so that there would be large fluctuations in this 
domain. We call this domain transition region hereafter. 

This circumstance is understood more clearly by study- 
ing the equation of state; the pressure is illustrated as a 
function of the inverse density as shown in Fig. 3 with 
T = 30MeV. At once it is noticed that this equation of 
state is similar to the van der Waals' one. This behav- 
ior is based on the fact that an interaction between two 



FIG. 4: The effective mass M as a function of the chemical 
potential at T = 30MeV. 



particles in many-particle systems is attractive at large 
distance and repulsive at short distance. In the present 
case, our system has both the attractive force (quark- 
quark interaction) and the repulsive one due to the Pauli 
principle. Hence our equation of state in the Fig. 3 is nat- 
ural. The straight line AB where the chemical potential 
is constant is the Maxwell constraction. If the hadron 
phase is regarded as "liquid" and the quark gas phase as 
"gas" , the present phase transition is opposite to that in 
the usual liquid-gas transition. The state "A" (and out- 
side) corresponds to the hadron (liquid) phase and the 
state "B" (and inside) to the quark gas phase. The seg- 
ment AC (BD) of the curve represents metastable state, 
which is called superheated liquid (supercooled gas). The 
remaining segment CD denotes unstable state. Therefore 
the quark matter jumps from A to B across the critical 
temperature T c . This behavior is also understood by the 
solutions of the gap equation (8) as illustrated in Fig. 4. 
Here three solutions are drawn as a function of /1 where 
four points (A, B, C and D) correspond to those in Fig. 3. 

Moreover we show the equations of state at three typ- 
ical temperatures in Fig. 5. The equation of state at 
T = lOOMeV is monotonic since the phase transition is 
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FIG. 5: The equation of states at T = 100, 60, 30MeV. 
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FIG. 6: The specific heat of the quark matter as a function 
of the temperature at /i = 0.25, 0.29, 0.31GeV. 



FIG. 7: The compressibility of the quark matter as a function 
of the temperature at /_t — 0.25, 0.29, 0.31GeV. The values 
at ft, = 0.31GeV are divided by 10. 



V. DISCUSSIONS AND SUMMARY 

It has been shown that there appears the large en- 
hancement of the specific heat and the compressibility 
near the tricritical point. These behaviors are due to 
the large fluctuation of the order parameter. Here we 
will clarify them in the context of the Landau mean field 
theory of phase transitions. 

To this end we take a\ as an order parameter of the chi- 
ral phase transition and investigate the thermodynamic 
potential (6) as a function of <j\. The other mean field 
02 is regarded as a function of <j\ through Eq.(8). If <j\ is 
denoted by 77, the thermodynamic potential is expanded 
as a power series of 77 as follows: 



second order and the others are first order. The equation 
of state at T = has been already calculated in Ref . |22^ . 

Next let us calculate the specific heat and compress- 
ibility by difference approximation of Eqs.(12) and (14) 
respectively. They are shown in Figs. 6 and 7 as a function 
of T at three chemical potentials: The case of /j = 0.25 
corresponds to the phase transition of the second order 
and the others to that of the first order. 

The black circles in these figures represent the calcu- 
lated values in the broken phase and the white ones in the 
symmetric one. Generally the specific heat in the broken 
phase is larger than that in the symmetric phase as well 
as the compressibility. The gray circles denote those in 
the metastable states. Although these are not realized 
within the mean field approximation, they may play an 
important role in the actual process. In fact soft modes 
will grow up in this transition region 18]- 23]. Namely 
the thermodynamic quantities of second derivatives are 
divergent near the temperature Ti, in particular the tri- 
critical point. These properties may be used as a signal of 
the existence of the phase transition of the quark matter 




u = w + An 2 + Bif + Dr 1 6 + ---, (15) 

where A, B and D are functions of T and fi. Moreover 
the D is assumed to be always positive because of the 
stability of our system. 

The order parameter should be determined by the sta- 
tionary condition: (dej/drf) = 2ri(A + 2Bri 2 + 3Drj 4 ) = 0. 
This equation has three solutions, 

(Vo = 0, 

1 _2 - -B±y/B 2 -3AD \ W ) 
I V± = 3D ' 

where rj± > is assumed. Here our discussion is sepa- 
rated into two cases according to the sign of B: (a) B > 
and (b) B < 0. In the case (a), the transition is second 
order. The stationary point with the minimum energy 
becomes 77 = 770 (for A > 0) or rj = r/ + (for A < 0). Since 
the critical phenomena of the second order transition is 
well known, we will not consider this case hereafter. In 
later discussion we will restrict ourselves to the other case 
(b) that is the first order transition. 

Then we have the following solutions depending on A, 
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B and D or (T, fx): 

(1) (r; ), (for A < 0) 

(2) 770, r; + , (for < A < B 2 /3D) (17) 

(3) 770 (for A > B 2 /3D) 

The solution in the parentheses means an unstable so- 
lution. In the case (2) which is called the transition re- 
gion, one of two extrema gives a stable state and the 
other a metastable state. The condition of the degen- 
eracy (tu(i]o) — Ld{i]+) or A — B 2 /AD) determines the 
critical temperature T c . The T\ (T2) is determined by 
A = B 2 /3-D (A = 0). These circumstances are confirmed 
in the previous section by calculating our selfconsistency 
equations (8) directly. An example is shown in Fig.l 
(/i = 0.31GeV). It is seen that the above case (2) is 
realized in 26 < T < 47MeV. 

Next we discuss the critical behavior of the specific 
heat near the temperature T±. Since the condition of 
T — T\ is given by A — B 2 /3D, we can expand as A — 
(B 2 /3D) = a(T — Xi) where a is a positive function of /j. 
Then the thermodynamic potential becomes 

u = ^-7^ m {X + B)(X-B) 2 , (18) 



where X = \/ B 2 — 3AD. Differentiating with respect to 
T, we get 



S = S» + —{B-^3aD{T 1 -T)), 



C = Co + ^^T^-T)-^ 2 . (19) 

It is seen that the entropy has finite gap across the tem- 
perature T\. This generates the latent heat as discussed 
in Ref. 16]. On the other hand, the specific heat di- 
verges at Ti, which is consistent with the numerical result 
(Fig. 6). As for the compressibility, we pay attention to 
the equation of state drawn in Fig. 3. Noting that the 
compressibility can be rewritten as k = —p(dp~ 1 /dP), 
it is related to the derivative of the equation of state. 
From the curve we see the compressibility is enhanced in 
the metastable phase, in particular divergent at T = T\, 
which is also seen in Fig. 7. 

The divergence of the specific heat and the compress- 
ibility can be also explained by using the thermodynamic 



potential. The temperature Ti is defined by B 2 — 3 AD, 
which gives 

u = luq + {A + Btf + )n\ + Dr)\ = lu + Drf + . (20) 

This cancellation of the second and fourth order terms 
leads to the large fluctuation of the order parameter. 

In conclusion, we have developed the chiral phase tran- 
sition of the first order in the framework of the NJL 
model. The quark matter discussed in this paper contains 
six kinds (two-flavor and three-color) of quark. Each of 
them has the same particle number, which means that 
our state is color singlet but not electric neutral. If this 
theory is applied to the bulk system, the charge neu- 
trality condition would be important. Our calculation 
has been carried out using the mean field approximation, 
which does not take into account the pairing correlation, 
color superconductivity [^-j^]. This effect may be im- 
portant for the high density and low temperature region 
and one of the future problems [^-jU]- We have ob- 
tained the equation of state of the quark matter, which 
is similar to the van der Waals' equation. Moreover the 
specific heat and the compressibility have been calcu- 
lated. They are enhanced in the broken phase and are 
divergent in the metastable state, in particular at the 
tricritical point. It may be expected that this singular- 
behavior appears in the expanding (cooling) process of 
the QGP generated in the high energy heavy-ion colli- 
sions. This expansion is like to the Big Bang and the in- 
flationary expansion of the early universe |32|-|34j. This 
singular behavior near the tricritical point may become 
a precursor of the chiral phase transition. 
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